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ABSTRACT 

A significant fraction of nearby elliptical galaxies are known to have high-density gas disks in their 
circumnuclear (CN) region (0.1 to a few kpc). Yet, ellipticals, especially luminous ones, show little 
signs of recent star formation. To investigate the possible cause of the dearth of star formation in 
these systems, we study the gravitational stability of CN gas disks embedded within the gravitational 
potentials of both the stellar bulge component and the central massive black hole (BH) in elliptical 
galaxies. We find that CN disks in higher mass galaxies are generally more stable than those in 
lower mass galaxies, because higher mass galaxies tend to have more massive BHs and more centrally 
concentrated stellar density profiles. We also consider the case in which the central stellar density 
profile has a core, which is often observed for ellipticals whose total stellar mass is higher than about 
10"'^"'^ Mq. Such a cored stellar density profile leads to more unstable CN disks than the power-law 
stellar density profile characteristic of less massive galaxies. However, the more massive BHs in high- 
mass galaxies act to stabilize the CN disk. Our results demonstrate that the gravitational potentials 
of both the central BH and the stellar component should be taken into account when studying the 
properties of CN disks, as their stability is sensitive to both the BH mass and the stellar density 
profile. Our results could explain the observed trend that less luminous elliptical galaxies have a 
greater tendency to exhibit ongoing star formation than giant ellipticals. 

Subject headings: galaxies: kinematics and dynamics — galaxies: formation — galaxies: stellar content 



1. INTRODUCTION 

Elliptical galaxies were believed to be completely dor- 
mant and therefore to consist of old stellar populations, 
due to the galactic wind which had blown the gas out 
and stopped star formation (Johnson & Axford 1971; 
Mathews & Baker 1971; Faber & Gallagher 1976). The 
discovery of an X-ray-emitting hot interstellar medium 
(ISM) in these objects has dramatically renewed our view 
of their physical properties and formation history (e.g., 
Forman et al. 1979, 1985). Warm gas is also detected in 
a significant fraction of ellipticals (55%— 60%), although 
the estimated mass is small (10'^ — lO'' Mq) (Phillips et al. 
1986; Ho et al. 1997a). In addition, multiwavelength ob- 
servations from optical to radio reveal the existence of 
the fair amount of H I gas (e.g., Knapp et al. 1985; Breg- 
man et al. 1992; Huchtmeier 1994; Knapp 1999; Morganti 
et al. 2006), dust (e.g., Knapp et al. 1989; van Dokkum & 
Franx 1995; Wiklind & Henkel 1995; Tomita et al. 2000; 
Colbert et al. 2001; Tran et al. 2001; Martel et al. 2004), 
and molecular gas (e.g., Wikhnd et al. 1995; Knapp & 
Rupen 1996; Vila-Vilaro et al. 2003; Young 2005; Nakan- 
ishi et al. 2007) in a significant fraction of ellipticals. Par- 
ticularly, some of the cold gas appears as a disk at the 
center of the galaxies. High-resolution images with Hub- 
ble Space Telescope (HST) uncovered that the majority 
of elliptical galaxies possess such dusty cold disks, with 
a typical mass of 10"* — 10^ Mq and sizes of 100 — 500 pc 
(e.g., Tomita et al. 2000). Molecular gas disks are also 
observed in the central region, with estimated masses of 
10^ — 10^ Mq and extents up to a few kpc (Knapp 1999; 
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Okuda et al. 2005; Young 2005; Nakanishi et al. 2007). In 
this paper we call such observed central cold structures 
"circumnuclear" (CN) disks. So far, there is no evidence 
of a correlation between the mass of the CN disk and 
the total stellar mass of ellipticals (e.g., van Dokkum & 
Franx 1995). 

The estimated densities of the CN disks are relatively 
high. CN disks in gas-rich disk galaxies obviously har- 
bor star formation (e.g., Phillips et al. 1986; Ho et al. 
1997b; Koda et al. 2005). However, it is still a matter 
of debate whether or not the CN disks in ellipticals have 
star formation. Ho et al. (1997a, 2003) and Ho (2005) 
suggest that there is no clear evidence of star formation 
in the central region of bright elliptical galaxies. On the 
other hand, Phillips et al. (1986) show that less luminous 
ellipticals tend to have star formation. It also seems that 
the properties of the ISM of ellipticals depend on galaxy 
mass: the detection rates of H I and molecular gas are 
higher in less luminous galaxies (Lake & Schommer 1984; 
Lees et al. 1991; Sadler et al. 2000; Sage et al. 2007). 
Okuda et al. (2005) measured the density and rotation of 
the CN disk in the radio galaxy 3C 31, and suggest that 
its disk is stable against local gravitational instability, 
consistent with the lack of evidence for star formation 
in this system (Owen et al. 1990). Young (2005) also 
performed a similar analysis for NGC 83 and NGC 2320. 
Therefore, the gravitational stability of the disk may be 
a key factor for the star formation condition in the CN 
disk (Kennicutt 1989) (but see also Koda et al. 2005). 

In the central region of elliptical galaxies, there are 
two empirical trends as a function of their mass. One 
is a well-established correlation between the black hole 
(BH) mass and the stellar mass of bulge (Magorrian et al. 
1998; Marconi & Hunt 2003; Haring & Rix 2004). The 
other one is the fact that the central stellar density pro- 
file loosely depends on the mass of galaxies (Faber et al. 
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1997; Ravindranath ct al. 2001; Stiavelli et al. 2001; Fer- 
rarese et al. 2006; Laucr ct al. 2007). The high-resolution 
images from the ffST enable the measurement of the stel- 
lar density profile in the central region of galaxies (Lauer 
et al. 1991, 1992b,a, 1993; Crane et al. 1993; Grillmair 
et al. 1994; Jaffe et al. 1994; Lauer et al. 1995). Lauer 
ct al. (1995) suggest two different categories for the in- 
ner profile of elliptical galaxies: luminous galaxies have 
a "core" and less luminous galaxies have a "power-law" 
profile without any evidence of a core. On the other 
hand, Ravindranath et al. (2001) and Rest et al. (2001) 
found galaxies with "intermediate" inner profiles that 
cannot be categorized as either cores or power laws. Tru- 
jillo et al. (2004b) attempt to link the central and global 
profiles of early-type galaxies by noting that the outer 
profiles of ellipticals are well known to be fitted by the 
Sersic law (r^/" law Sersic 1968), and whose index n 
is well-correlated with the global properties of galaxies, 
such as the effective radius and the total luminosity (e.g., 
Caon et al. 1993; Graham & Guzman 2003). They sug- 
gest that the inner power-law profile can be explained 
as the extension of the outer Sersic profile, except for 
galaxies with central cores. Therefore, among power-law 
galaxies, the brighter members have profiles with larger 
n, and hence steeper power-law slope for the inner pro- 
file. On the other hand, Trujillo et al. (2004b) propose 
that luminous galaxies with cores can be described with 
a "core-Sersic" profile, which is a combination of a core, 
approximated as an inner power-law profile with a shal- 
low slope, and an outer Sersic profile. 

Motivated by these facts, we study the gravitational 
stability of the CN disk for elliptical galaxies with differ- 
ent masses. The stability of a rotating disk can be de- 
scribed by the so-called Toomrc's Q-parameter (Toomre 
1964; Goldreich & Lynden-Bell 1965). This parameter 
also provides a criterion for star formation in rotating 
disks (e.g., Kcnnicutt 1989). Tan & Blackman (2005) 
study the stability of the CN disk analytically. They con- 
clude that CN disks are inevitably unstable, and should 
have star formation activity even in giant ellipticals, such 
as M87. However, so far, no indication of star formation 
has been observed in M87. Although Tan et al. (2007) 
report the detection of the molecular gas in the central 
region of M87, the optical spectral features of M87's CN 
disk are not consistent with star formation. The nu- 
clear spectrum of M87 shows classical signatures of low- 
ionization nuclear emission- line regions (Ho et al. 1997c), 
which as a group are consistent with being active galac- 
tic nuclei (AGNs) accreting at a low rate (Ho et al. 2003; 
Ho 2004). 

The gravitational potential in the CN region is gov- 
erned by the central BH as well as the stars in the inner 
bulge. The study of Tan & Blackman (2005), however, 
ignores the potential of the stellar bulge component. Our 
analysis of the stability of the CN disk in ellipticals takes 
into account both contributions to the potential and their 
dependence on the global mass of the system. Section 2 
summarizes our method to analyze the stability and de- 
scribes models of the CN disk and the stellar potential 
in the CN region. Section 3 shows the results, and in 
Section 4 we discuss how our results fit with the recent 
observed properties of the central region of ellipticals. 

2. METHOD 



The stability of a self-gravitating disk can be analyzed 
by Toomrc's Q-paramctcr (Toomre 1964; Goldreich & 
Lynden-Bell 1965). The definition of Toomre's Q pa- 
rameter is 



where Cg, k, and arc the sound velocity, epicycle fre- 
quency, and surface density of the gas disk. If Q < 1) 
the disk is unstable. The specific frequency for a disk is 
described as 

n' = R^ + 4n\ (2) 

where Q is the circular frequency. 

We assume that the disk is steady, i.e. the accretion 
rate is the same at different radii. The accretion rate can 

be written as 

M{R) = ~2TTRa^ccCsMRY /n{Rf]T,diR) = const., 

(3) 

where aacc is the dimensionless viscosity parameter 
(Shakura & Smiyaev 1973; Pringle 1981; Frank et al. 
2002). This allows us to derive the density profile of 
the CN disk, Sd(i?), once the total mass, Md, and the 
radius, Rd, of the disk are fixed. The circular density 
profile ri(i?) is determined by the gravitational poten- 
tial. For simplicity, we assume that the gravitational 
potential is dominated by the central BH and the stel- 
lar component, and that the contribution from the CN 
disk is negligible. As will be shown in Figure 1 below, 
this assumption is valid, except for the cases of massive 
and compact CN disk in smaller mass galaxies, which 
are inevitably unstable. We also assume that the stellar 
density is much higher than the dark matter density in 
the central region. 

We examine the stability of the CN disk within the 
potential of the central region of spherical galaxies with 
different total masses of the stellar bulge, Mg. We adopt 
the relation between the BH mass and the stellar mass 
from Hiiring & Rix (2004), 

log(MBH/M0) = 8.20 + 1.121og(M,/lO"M0). (4) 

Hence, once Mg is fixed, we can calculate the potential 
from the BH with the mass of Mbh- 

As mentioned in Section 1, the mass dependence of 
the inner stellar density profile is discussed by a number 
of authors. The first thorough study was carried out 
by the "Nukcr team," who introduced the "Nuker law" 
to describe the observed surface brightness profile of the 
central regions of galaxies. The Nuker law (Lauer et al. 
1995; Faber et al. 1997) is given by 

I{R) = 7b2^''-'^^/"(i?/i?b)"''[l + (i?/i?b)"]^'^"''^/". (5) 

The asymptotic logarithmic slope inside i?b is 7; 
the asymptotic logarithmic outer slope is f3; and a 
parametrizes the sharpness of the break. The break ra- 
dius, i?b, is the point of maximum curvature in log-log 
coordinates. The break surface brightness, /b, is the sur- 
face brightness at i?b- This function is designed to fit the 
surface brightness profile in the inner region, and not to 
describe the entire profile. Graham et al. (2003) pro- 
posed a new formula that is a combination of an inner 
power-law profile and an outer Sersic law. They call it 
the "core-Sersic law," which is described as 

IiR) = I' [1 + (i?b/ii)T^" 
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with 



exp{-M(i?« + i?g)/i?«]^/("=«)}, (6) 

/' = /b2-^/" exp [62V("=")(i?b/i?e)^/"^] • (7) 

The parameter R^, is the break radius to separate the 
inner power law with a slope of 7 from the outer Sersic 
law with effective radius i?e and index rig, and /b is the 
surface brightness at Rh- The parameter a controls the 
sharpness of the transition between the inner and outer 
profiles, where a higher value leads to sharper transitions. 
Trujillo et al. (2004b) suggest that the core-Sersic law 
can be simplified with a — » 00 and still provide a good 
description to the observed profiles (see also Ferrarese 
et al. 2006). In this limit, 

I{R) = h[{Rb/Ryu{Rb- R) 

where u{x — a) is the Heavisidc step function. 

It is still controversial which of the two formalisms bet- 
ter describes the surface brightness profile in the cen- 
tral region of galaxies (Ferrarese et al. 2006; Lauer et al. 
2007). In this paper, we adopt the corc-Scrsic law simply 
for computational convenience to link the density profiles 
in the inner region to the total mass of the stellar bulge. 
Although the simplified version of the core-Sersic law as 
given in equation (8) has been adopted to fit recent ob- 
servations (Trujillo et al. 2004b; Ferrarese et al. 2006), 
wc use equation (6) by setting a to 5.0, for computa- 
tional convenience. Lauer et al. (2007) demonstrate how 
the sharp transition generated from equation (8) leads 
to poor fits to the observed surface brightness profiles. 
The adopted a provides a less sharp transition. This 
assumption also guarantees that the deprojected density 
profile does not increase with radius with our adopted 
parameters (see below), but the profiles is still close to 
the simplified formula of equation (8) used in the obser- 
vations. 

As a comparison between the Nuker law (eq. 5) and 
the core-Sersic law (eq. 6), we briefiy mention how the 
logarithmic gradient, j'{R'), at radius of R' can be de- 
scribed for each fitting function (see also Trujillo et al. 
2004b). Rest et al. (2001) show that for the Nuker law 



7'(i?') 



dlogI{R) 



^ + f3{R'/Rk) 

1 + {R'/R^Y 



dlogR 

For the core-Sersic law (eq. 6), Trujillo et al. 
show that the slope can be written as 



-. (9) 
(2004b) 



where /tot is the total bulge luminosity of the integration 
for < .T < 00 (sec also Graham & Driver 2005, and 
references therein). For simplicity, we define the total 
bulge mass, Mg, as the total mass of the integration of the 
Sersic law of equation (11), regardless of the existence of 
the core. Then, assuming a constant mass-to-luminosity 
ratio T at the different radii, we can describe the surface 
brightness at the break radius, /b, for the core-Sersic law 
(eq. 6) as 

]\/f J,2na 

7j^ = i^__%__e-''4iib/iJe)^/-. (14) 



j'(R') = — (l/i?e)^/"=i?'"(i?'"+i?JJ) 
7ls 



l/Ka)-l + ^CRbW 



1 + {Rb/R'Y 
(10) 

The core-Sersic law (eq. 6) is equivalent to the Sersic 
law, when i?b = and 7 = 0. Then, the Sersic law is 
described as 

/(x)=AsCxp(-M'/"=), (11) 
with X = R/Rg. Prugniel & Simien (1997) derive the 
numerical solutions 

-4. = /tot . o (12) 



bs = 2ns-\+ 0.009876/ns 



(13) 



T 27rnsr(2ns) 
It is well known that Re and rig correlate with the to- 
tal luminosity of galaxies (e.g., Kormendy 1977: Caon 
et al. 1993; Graham 2002). From their fits of T/ST im- 
ages of Virgo early- type galaxies using equation (8), Fer- 
rarese et al. (2006) find that the profiles of bright galaxies 
{Mb < —20.5 mag) require a core, while those of less lu- 
minous galaxies can be described with the Sersic law (see 
also Graham & Guzman 2003). They also derive the fol- 
lowing relations between the g-band effective radius and 
B-band absolute magnitude: 

logi?e = -0.055(AfB + 18) + 1.14 (15) 
for the Sersic galaxies, and 

logi?c = -0.22(AfB-|- 18) + 1.5 (16) 
for the core-Sersic galaxies (Ferrarese et al. 2006, for error 
estimates, see). The units of _Rc is arcseconds, which for 
Virgo corresponds to 1" = 80.1 pc. In addition, they 
suggest that the Sersic index ris and the B-band absolute 
magnitude are related by 

log Us = -0.10(Mb + 18) + 0.39, (17) 
valid for galaxies of both classes. In this paper we adopt 
these relations for computational convenience. Note that 
these relations are still not well established. For exam- 
ple, Graham & Guzman (2003) and Graham et al. (2006) 
suggest a curved relation between the luminosity and ef- 
fective radius rather than two power laws. 

We further assume that the stellar mass-to-light ra- 
tio for early-type galaxies in the B-band is Tb = 
Ms/ {Lb / Lb,q) = 7, with Mb.q = 5.48 mag, and that 
Tb is independent of mass; these assumptions are suf- 
ficiently accurate for early-type galaxies (Trujillo et al. 
2004a). With some exceptions, Ferrarese et al. (2006) 
find that there is a critical luminosity that separates 
the Sersic galaxies from the core-Sersic galaxies. With 
the assumptions made above, their critical luminosity of 
Mb,c = -20.5 mag corresponds to Ms^c = 1-73 x 10^^ 
Mq; systems with masses lower than Ms,c are assumed to 
follow the Sersic law. Then, we derive the half-mass ra- 
dius using equations (15) and (16), with the assmnption 
■that the g-hand half-light radius is similar to the half- 
mass radius. Once we fix Mg, all the parameters of equa- 
tion (6) can be determined, and equation (6) provides 
the projected mass density profile of Yss{R) = Ib{R)Tb- 
Assuming spherical symmetry, we can derive the three- 
dimensional stellar mass density Ps{r) through (e.g., Bin- 
ney & Tremaine 1987) 

1 r dMR) dR 

The analytical formula of three-dimensional stellar den- 
sity profiles is also discussed in Prugniel & Simien (1997), 
Mazure & Capelato (2002), Terzic & Graham (2005), and 
references therein. 
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Fig. 1. — Stellar density profile as a function of three-dimensional 
radius (top) and the total mass of stars and BH within the radius 
{bottom) for four different models indicated in the upper panel. 



TABLE 1 

Stellar Bulge Density Profiles for the Host Galaxies 



Model 




Profik- 


Bo 




By, 




Mb"" 




(M . ) 




(kpc) 




UH-) 




(mag) 


SIO 


10"' 


Scrsic 


1.03 


2.14 






-17.4 


Sll 


lOii 


Sersic 


1.41 


3.81 






-19.9 


cSll 


10" 


core-Sersic 


6.66 


3.81 


9.32 


0.1 


-19.9 


cS12 


1012 


core-Sersic 


23.6 


6.77 


331 


0.1 


-22.4 



''B-band luminosity estimated with a stellar mass-to-light ratio of 
Ts = M^/{Lb/Lb,q) = 7. 



3. RESULTS 

3.1. Stability of CN Disks with Fixed Mass and Radius 

As mentioned in Section 1, so far no correlation has 
been found between the mass or size of CN disks and 
the host galaxy mass. In this section we study how the 
stability of CN disks with fixed mass and radius depends 
on the total mass of the host galaxies. Here, we focus on 
a disk with mass log(Afd/M0) = 8 and radius i?d = 2.5 
kpc. which is roughly the same as the CN disk observed 
in NGC 4476 by Young (2002). Disks with a range of 
difi^erent masses and radii will be discussed in the next 
section. Note that once Md and Rd are fixed, Sd(^) 
does not depend on aacc or Cg. Hence, Q{R) described 
by equation (1) is independent of aacc and is simply pro- 
portional to Cs = i/7dfcBT'd/(A«mp), where 7d, fce, T^, 
ji, and TTip are the specific heat, Boltzmann's constant, 
the gas disk temperature, the mean molecular mass, and 
the proton mass, respectively. Throughout the paper, we 
fixed Cs to a value corresponding to /U = 0.6, 7d = 1, and 
Td = 30 K, which is a typical temperature of observed 
CN disks (e.g., Wiklind et al. 1995; Tran et al. 2001). 

Stability of the CN disk is examined for host galaxy 
models with stellar bulge masses of log(Ms/M0) = 10, 



11, and 12. We apply the Sersic law to the model with 
log(Afs/Af0) = 10 (model SIO), and the core-Sersic law 
to the model with log(Ms/MQ) = 12 (model cS12). Since 
the adopted scaling relation used in Section 2 is known to 
have significant scatter and the model of log(Ms/M0) = 
11 is close to the critical mass of Mg^c = 1-73 x 10^^ Mq 
for the two stellar density profiles, we consider two cases 
of the Sersic (model Sll) law and the core-Sersic (model 
cSll) laws for host galaxies with log(Ms/M0) = 11. All 
the parameter values for the stellar density profiles and 
the estimated B-h&nd absolute magnitude are summa- 
rized in Table 1. We choose 7 = 0.1 for the core-Sersic 
law of equation (6), which is a relatively low value among 
observed cores (e.g., Ferrarese et al. 2006); this is a con- 
servative choice because it leads to a maximally unstable 
CN disk. 

Figure 1 shows that the stellar density profile and the 
total mass of stars and BH as a function of radius for the 
four models. Comparison between models Sll and SIO 
shows that the density profile is steeper for higher mass 
galaxies, which have a larger ng. The core-Sersic models 
lead to a lower density in the inner region. Consequently, 
model cS12 has even lower density than model Sll and 
a similar density to model SIO within the break radius of 
cS12. However, higher mass galaxies have a more massive 
BH. As a result, the total mass within a given radius 
is generally higher for the higher mass galaxies, except 
around the break radius for the core-Sersic models. 

Figure 2 presents the values of Toomre's Q-parameter 
at difi^erent radii of the CN disk for the four models. 
We also examine two additional cases for each model; 
one is the case assuming no BH (dotted line), and the 
other one (dashed line) ignores the stellar potential (i.e., 
only the; Keplerian potential of the BH is considered). 
As described in Section 2, the circular velocity (Fig. 3) 
is calculated purely by the BH and/or stellar potential, 
and the gas disk density profile (Fig. 4) is derived from 
the condition of M{R) — constant in equation (3). 

The BH-only potential results (dashed lines) demon- 
strate that since higher mass galaxies have a more mas- 
sive BH, CN disks are more stable in higher mass galax- 
ies. However, for the assumed CN disk, even model cS12 
still has Q < 1 — i.e., the CN disk is unstable. The 
stellar -|-BH potential results clearly demonstrate that 
the stellar potential greatly helps to stabilize the CN 
disk. Moreover, comparison between the results for the 
stellar+BH potential and the results of the case of stellar- 
only potential indicates that the BH helps to stabilize the 
CN disk only in the inner region (a few 100 pc, depending 
on the BH mass and the stellar density profiles). There- 
fore, the stellar potential cannot be ignored in studying 
the dynamics of the CN disk. 

For the Sersic models, the stellar potential stabilizes 
the CN disk more in the higher mass galaxies, because 
higher mass galaxies have larger rig, i.e. higher central 
concentration of the stellar potential (see also Fig. 1). At 
a fixed stellar mass of log(Ms/M0) = 11 the core-Sersic 
model leads to a lower Q- value and a more unstable disk. 
The surface density profile of the CN disk is sensitive to 
this change of the stellar potential profile (Fig. 4). The 
stellar-only potential results clearly demonstrate that the 
disk gas density has a peak around the break radius. Due 
to this peak density, the Q-value becomes as low as the 
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Fig. 2. — Toomrc's Q-value as a function of radius for the CN disk 
with M(j = 10* Mq and = 2.5 kpc in galaxy models indicated 
in the panels. The black solid line presents the Q-value calculated 
by taking into account both the BH and the stellar potential. The 
dotted line only includes the stellar potential, while the dashed line 
takes into account only the potential of the central Bll, which is 
correlated with the mass of the stellar component as assumed in 
equation (4). The gray solid line marks Q = 1. 

case of the BH-only potential. However, the BH can 
stabilize the disk in the inner region, and the stellar -|-BH 
potential case is much more stable. As a result, model 
cSll is more stable than model SIO. 

Note that, as shown in Lauer et al. (2007), the Sersic 
law often underestimates the surface brightness in the 
central region where a high-density stellar component is 
often seen. Ferrarese et al. (2006) invoke an additional 
component, which they call "nuclei," to fit this central 
high-density component. Although we do not include the 
nuclei in this study for simplicity, such a compact central 
component also stabilizes the CN disk in a similar way to 
the BH. Cote et al. (2006) find that the estimated mass 
of the nuclei is similar to what the BH mass and bulge 
mass relation, such as equation (4), predicts. Therefore, 
the effect of the nuclei would be similar to assuming a 
factor of two more massive BH. 

In general, CN disks in higher mass galaxies are more 
stabilized, due to higher mass BH and more centrally 
concentrated profile of the stellar component. However, 
if the host galaxy has the core-Scrsic law profile, the CN 
disk is more unstable, compared with the CN disk in 
a galaxy with the same stellar mass and BH mass, but 
having Sersic profile. 

3.2. Stability of CN Disks with Various Masses and 

Radii 

In this section, we again consider the four models 
shown in Table 1, but study the stability of the CN disks 
with different masses and radii. Figure 5 shows the Q- 
value for the CN disks with radii in the range of 0.05 — 15 
kpc and masses spanning log(M(i/M0) =4—10. Here, 
the figure shows the minimum Q-value in each model. 
For example, Qmin ~ 3.15 at ra = 380 pc for model cS12 
in Figure 2. 

As expected, CN disks with higher gas mass and 
smaller radii are more unstable. Comparison between 
SIO and Sll demonstrates that for two CN disks with 
the same mass and radius in a Sersic-law galaxy, the disk 
in the higher mass galaxy is always more stable, due to 
the higher mass BH and the more centrally concentrated 
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Fig. 3.-- The circular velocity, t)c, as a function of the radius of 
the CN disk with = 10* Mq and = 2.5 kpc in galaxy models 
indicated in the panels. The solid line takes into account both the 
BH and the stellar potential. The dotted line only includes the 
stellar potential, while the dashed line takes into account only the 
potential of the central BH. 
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Fig. 4. — The surface density profile, Sd(-R), of the CN disk with 
M(j = 10* Mq and R^ = 2.5 kpc in galaxy models indicated in 
the panels. The solid line takes into account both the BH and the 
stellar potential. The dotted line only includes the stellar potential, 
while the dashed line takes into account only the potential of the 
central BH. 

stellar density profile. If the CN disk is smaller than the 
break radius, the highest mass model (i.e. model cS12) 
is the most stable. For a CN disk with a size compara- 
ble to the break radius, model cS12 leads to a less stable 
condition than model Sll, and for the larger CN disk 
model cS12 is as stable as model Sll. Therefore, for 
large CN disks in high mass galaxies, the size of the disk 
with respect to the break radius of the host galaxies is 
an important factor for the stability. 

4. DISCUSSION AND CONCLUSIONS 

Our analysis of the stability of CN disks reveals the 
following trends. 

(1) The CN disk is stabilized by the presence of the 
central BH. Since higher mass galaxies tend to have 
more massive BHs, the CN disk is more stable in 
higher mass galaxies in general. 

(2) The stellar potential is also important for the sta- 
bility of the CN disk. 
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Fig. 5. — The minimum Q-valuc for the CN disk with diflfcr- 
ent masses and radii. The upper-left, upper-right, lower-left, and 
lower-right panels show the results of model SIO, Sll, cSll, and 
cS12, respectively. The thick solid line indicates (5 = 1. The dot- 
ted lines and thin solid lines correspond to contours for Q < 1 and 
Q > 1. The levels arc separated by ilogQ = 1.0. Open circles 
denote CN disks observed in CO emission by Young (2002, 2005). 
The rest of the symbols show the mass and size of the CN disk es- 
timated from the optical color excess of dusty disks observed with 
HST (Tomita ct al. 2000). Solid triangles and open squares cor- 
respond to the galaxies whose profile has been fitted with a Sersic 
law and a core-Sersic law, respectively (Ferrarese et al. 2006; Lauer 
et al. 2007); open diamonds mark galaxies with unclassified profiles. 
NGC 4476, which has been observed both in CO (Young 2002) and 
with HST (Tomita et al. 2000), is connected with a dotted line. 

(3) For Sersic-law galaxies, the CN disk is more stable 
in higher mass galaxies because they tend to have 
more centrally concentrated stellar density profiles. 

(4) The existence of a central stellar core in luminous 
ellipticals makes the CN disk unstable, especially 
around the break radius. 

As discussed in Section 2, the central surface bright- 
ness profiles for relatively low-mass galaxies (Mg < 10"^^ 
Mq) arc generally described by the Sersic law, with the 
trend of higher mass galaxies tending to have larger ng. 
Points (l)-(3) above indicate that CN disks in lower mass 
galaxies are more unstable. On the other hand, rela- 
tively high-mass galaxies (log(Afs/M0) > 11) tend to 
have stellar density profiles that contain a central core, 
which leads to a more unstable CN disk compared to the 
Sersic law. On the other hand, since such systems also 
have more massive BHs, CN disks in higher mass galax- 
ies are kept stable, especially within the break radius. 
These trends can explain the observational trends out- 
lined in Section 1. Because CN disks are more stable in 
more luminous ellipticals, this provides a natural expla- 
nation for central star formation to be curtailed in giant 
ellipticals, whereas less luminous ellipticals apparently 
have host nuclear star formation with greater case, even 
though cold ISM in the form of dusty nuclear disks are 
observed to be just as prevalent in both environments. 

It is also worth stressing that, as seen in Figure 4, the 
CN disk can remain stable even if the density of the disk 
is more than 100 Mq pc~^ within 100 pc. This is because 
the BH stabilizes the CN disk in the central region, as 
seen in Figure 2. This density is much higher than the 



canonical density threshold for star formation suggested 
by both observations (Kennicutt 1989) and theory (e.g., 
Schaye 2004), Sth ~ 3 - 10 Mq pc'^. Our stability 
analysis offers a simple explanation for the lack of star 
formation (e.g., Okuda et al. 2005; Young 2005) in CN 
disks that otherwise have high densities. 

At a fixed total stellar mass for the host galaxy, galax- 
ies with higher ng Sersic law have more stable CN disks, 
while the core-Sersic law leads to more unstable CN disks 
than the Sersic law. Therefore, the frequency of central 
star formation activity should depend on both the cen- 
tral stellar density profile as well as the BH mass. Re- 
cent near-ultraviolet observations performed with Galaxy 
Evolution Explorer (GALEX) indicate that some frac- 
tion of early-type galaxies have a small amount of recent 
star formation (Yi et al. 2005). Based on semi-analytic 
model predictions Schawinski ct al. (2006) demonstrate 
that the observed near-ultraviolet color distributions of 
early-type galaxies can be explained if there is a criti- 
cal BH mass at the fixed velocity dispersion for galaxies 
that have recent star formation. They suggested that the 
critical BH mass likely comes from strong AGN heating 
created by the massive BH (e.g., Binncy & Tabor 1995; 
Silk & Rees 1998; Kawata & Gibson 2005; Springel et al. 
2005; Croton et al. 2006; Bower et al. 2006; Ciotti & Os- 
triker 2007). However, our study demonstrates that the 
existence of a critical condition for star formation at the 
fixed velocity dispersion can be due to the stability of the 
CN disk, which is governed by both the stellar density 
profile and the BH mass. 

As mentioned in Section 1 , the size and mass of the CN 
disk have been measured by several observational stud- 
ies. Some of these measurements are plotted in Figure 
5 to compare with our model predictions. The circles in 
Figure 5 correspond to objects whose CN disk mass and 
radius were estimated from CO emission observed with 
interferometers in Young (2002, 2005). Figure 5 also con- 
tains objects whose disk properties were estimated from 
the color excess of dust features measured in optical im- 
ages of the central region of early-type galaxies observed 
with the HST (Tomita et al. 2000)*. The stellar density 
profiles of the central region for the majority of the galax- 
ies in Tomita ct al. (2000) have been studied by Ferrarese 
et al. (2006) and/or Lauer et al. (2007). This allows us 
to further distinguish the objects by their central profile 
type. 

Following Section 2, we estimated the stellar mass 
from the S-band luminosity, assuming a stellar mass- 
to-light ratio of = Mg/ [Lb / Lb,q) = 7. The S-band 
luminosities come from the NASA/IPAC Extragalactic 
Database (NED). We plot the data for galaxies with 
masses \og{Ms/MQ) < 10.5 and log(Ms/M0) > 11.5 in 
the upper-left and lower-right panels of the figure, re- 
spectively; the rest of galaxies are shown in both the 
upper-right and lower-left panels. 

Only one object — NGC 4476 — has been measured 
with both techniques (CO from Young (2002) and dust 
from Tomita et al. (2000)); the two data points are con- 
nected with the dotted line. Young (2002) obtained a 
much higher gas mass of 1.1 x 10^ Mq within 2.4 kpc, 

■* Tomita ct al. (2000) find that some galaxies in their sample 
have irregular dust morphology. Figure 5 only shows data for galax- 
ies whose dust morphology is classified a disk. 



Stability of Circumnuclear Disk in Ellipticals 



7 



compared with a mass of 1.0 X 10"^ Mq within 1.8 kpc 
in Tomita ct al. (2000). The order-of-magnitudc discrep- 
ancy in mass cannot be explained by the slightly different 
distances adopted by these authors (18 Mpc by Young 
(2002) and 16.8 Mpc by Tomita et al. (2000)). It is note- 
worthy that all of the CO-measured CN disks have sys- 
tematically larger disk masses than the dust-mcasurod 
systems, suggesting that at least part of the discrepancy 
may be due to systematic errors in the estimated gas 
masses. Both sets of observations require a conversion 
factor to arrive at a gas mass. In the case of the CO 
observations, a standard Galactic CO-to-H2 conversion 
factor was adopted. The CO-to-H2 conversion factor, 
however, may be systematically lower in regions of high 
metallicity (Arimoto et al. 1996), as is the case in the cen- 
tral regions of massive galaxies, a suggestion supported 
by radiative transfer calculations in numerical simula- 
tions of CN disks (Wada & Tomisaka 2005). If this is 
the case, then the gas masses from Young (2002, 2005) 
are overestimated. 

A similar caveat applies to the gas masses estimated 
from the dust extinction, which assume a Galactic gas- 
to-dust mass ratio (Tomita et al. 2000). However, dust 
mass estimated from optical extinction should be consid- 
ered as lower limit to the true dust masses (Goudfrooij 
& de Jong 1995; Tran et al. 2001; Martel et al. 2004). 
For example, Goudfrooij & dc Jong (1995) show that 
the dust masses derived from the far-infrared emission 
are roughly an order of magnitude higher than those es- 
timated from optical extinction alone. They argue that 
the optical extinction may underestimate the dust mass 
by about factor of 2 due to the assumption that the dust 
is in front of the stars (see also Martel et al. 2000, who 
claim that this effect is more significant). Goudfrooij & 
de Jong (1995) also discuss that the discrepancy in the 
estimated dust mass is mainly due to the presence of dif- 
fusely distributed dust, which cannot be seen as optical 
extinction. We note, however, that such a diffuse com- 
ponent is unlikely to be present as a rotating disk, as it 
is more likely to be supported by velocity dispersion. In 
any case, it seems plausible that the gas masses derived 
from the optical dust features may be an underestimate 
of the true gas mass, thereby narrowing the disagree- 
ment with the CO-based masses. Until more accurate 
gas masses are available, it is difficult to draw more quan- 
titative comparisons between our model predictions and 
observations. 

To sec if there is any sign of ongoing star formation in 
the sample of galaxies shown in Figure 5, we have done 
a careful search of the literature to inspect published op- 
tical spectra (Phillips et al. 1986; Bonatto et al. 1989; 
Ho ct al. 1997c; Cretton et al. 2000; Donzelh & Davoust 
2003). Galaxies with ongoing star formation show optical 
emission-line ratios that are readily distinguishable from 
other sources of ionization (e.g., active galactic nuclei: 



sec Ho ct al. 1997c). Among all the objects, only NGC 
4526 has a spectral classification consistent with stellar 
photoionization (Ho et al. 1997c). The rest either have 
no star formation (14 galaxies) or have insufficient spec- 
tral information to tell (5 galaxies). Although NGC 4526, 
which is highlighted in Figure 5, is in the stable (Q > 1) 
region, interestingly, NGC 4526 has one of the most un- 
stable CN disks among the iJS'T-measured sample. Also, 
none of the CO-measured CN disks, which are located in 
the unstable region in Figure 5, show a clear indication of 
star formation. Although this contradicts with our pre- 
diction, because of the above-mentioned ambiguities in 
the estimates of the gas masses from the observations, it 
is difficult to conclude if these galaxies require additional 
process to explain their suppression of star formation. 

It is clear from inspection of Figure 5 that once the 
mass of the CN disk becomes high enough, the disk will 
inevitably become unstable for star formation, one of the 
consequences of which may be to aid gas fueling to the 
AGN (Wada 2004). This type of situation is likely to 
be realized in the aftermath of a gas-rich major merger, 
whereby the large amount of the; gas dissipated toward 
the center can generate a high-mass CN disk. High- 
resolution millimeter observations of luminous infrared 
galaxies indeed suggest that systems with larger cen- 
tral gas surface densities tend to either form stars with 
greater efficiency or have a higher probability of hosting 
an AGN (Bryant & ScoviUe 1999). 

Finally, we speculate that star formation induced in an 
unstable CN disk may be directly related to the forma- 
tion of kinematically decoupled cores (KDCs) that are 
often seen in the centers of early- type galaxies (e.g., Ko- 
rmendy 1984; Bender 1988; Bender & Surma 1992; Caon 
et al. 2000) as well as the central stellar disk that many 
power-law ellipticals seem to have (Lauer et al. 2005). 
One of the intriguing clues regarding the origin of KDCs 
is that large cores arc always old (McDermid et al. 2006). 
If KDCs are a by-product of unstable CN disks, then 
large KDCs require large, massive CN disks, which prob- 
ably can only be formed through major, gas-rich mergers, 
which are most prevalent at earlier epochs. 
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